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1. Introduction 

Hydrodynamics is well a established effective approach to flow phenomena. Here we shall present its nonextensive 
relativistic version from the point of view of high energy collision physics [1-6] . The characteristic feature of such 
a processes is the production of a large number of secondaries (their multiplicities at present approach ~ 10 3 ). 
Already in 1953, when there were only ~ 10 particles produced and registered, mostly in cosmic rays experiments, 
it was found [1] that such processes can be very effectively described by using a thermodynamic approach in 
conjunction with hydrodynamical flow. In this so called Landau Hydrodynamic Model the secondaries secondaries 
were considered to be a product of decay of some kind of hadronic fluid produced in such collisions, which was 
expanding before hadronization [1]. This model has recently been updated to describe recent experimental data 
[2]. Since then there was a number of successful attempts to develop new solutions for both the Landau model [3] 



* E-mail: osada<Qph.ns.musashi-tech.ac.jp 
^ E-mail: wilk@fuw.edu.pl 



Noncxtensive perfect hydrodynamics - a model of dissipative relativistic hydrodynamics? 



and for the so called Hwa-Bjorken version of the hydrodynamical model [4], for which a new class of solutions has 
been found [5] . Hydrodynamic models of different types were therefore frequently used in the phenomenological 
description of multiparticle production processes at high energies, especially for high energy nuclear collisions 
[6]. These are of special interest to us (being currently investigated at RH1C at Brookhaven with the newly 
commissioned LHC at CERN joining soon) because it is widely believed that in such collisions a new form of 
hadronic matter, the so called Quark- Gluon- Plasma (QCD), will be produced [7]. 

So far all hydrodynamic models in this field were based on the usual Boltzmann-Gibbs (BG) form of statistical 
mechanics. The only works discussing some general features of the nonextensive hydrodynamics, which we are 
aware of [8], use a nonrelativistic approach and are therefore not suitable for the applications we are interested in. 
On the other hand it is known that an approach based on non-extensive statistical mechanics (used mainly in the 
form proposed by Tsallis [9] with only one new parameter, the nonextensivity parameter q) describes different 
sets of data in a better way than the usual statistical models based on BG statistics, cf., [9] for general examples 
and [10-f 5] for applications to multiparticle production processes. Roughly speaking, all observed effects amount 
to a broadening of the respective spectra of the observed secondaries (both in transverse momentum space 
and in the rapidity space), they take the form of g-exponents instead of the naively expected usual exponents: 
cxp(— X/T) exp {-X/T) = [1 - (1 - q)X/T] 1/(1 - q) . From these studies emerged a commonly accepted 
interpretation of the nonextensivity parameter q (in fact, \q — 1|) as the measure of some intrinsic fluctuations 
characteristic for the hadronizing systems under consideration [fO, 11]. For q > 1 and in the transverse 
momentum space it could be fluctuation of the temperature T corresponding to some specific heat parameter 
C. In this case q — 1 = C and therefore it should be inversely proportional to the volume of the interaction 
region. This effect is indeed observed [15]. In rapidity space these are fluctuations of the so called partition 
temperature, T v t = E/(n) [14], which are precisely the same fluctuations that lead to the Negative Binomial form 
of the observed multiplicity distributions, P(n; (n); k), with its characteristic parameter k given by k = l/(q — 1) 
[14, 16]. In the case of q < 1, the interpretation is not at present clear [11]. It seems that the first role of the 
parameter q is to restrict the allowed phase space [17]. Actually, the conjecture associating q with fluctuations 
has already been formalized as a new branch of statistical mechanics called superstatistics [18]. It should be 
noted that there are also arguments connecting nonextensivity with some special dynamical correlations existing 
in the system under consideration [19], but their connection with fluctuations is not yet fully clear (cf., [i6]) Gndnotcl 

The physical picture emerging from the above experience is that, instead of a strict local thermal equilibrium 
customarily assumed in all applications of statistical models (including hydrodynamic), one rather encounters 
a kind of stationary state, which already includes some interactions. It can be introduced in different ways. 
For example, in [20] it was a random distortion of the energy and momentum conservation caused by the 
surrounding system. This results in the emergence of some nonextensive equilibrium. In [21, 22] the two-body 
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energy composition is replaced by some generalized energy sum h(Ei,E2), which is assumed to be associative 
but which is not necessarily simple addition and contains contributions stemming from pair interaction (in 
the simplest case). It turns out that, under quite general assumptions about the function h, the division 
of the total energy among free particles can be done. Different forms of the function h lead then to differ- 
ent forms of the entropy formula, among which one encounters also the known Tsallis entropy. The origin 
of this kind of thinking can be traced back to the analysis of the q-Hagedorn model proposed some time ago in [23]. 

All phenomenological applications of hydrodynamic models to the recent multiparticle production processes 
show that, although perfect (nonviscous) hydrodynamics successfully describes most RHIC data [3, 24], there are 
indications that a hadronic fluid cannot be totally ideal. For example, the perfect fluid dynamical calculations 
with color glass condensate initial state could not reproduce the elliptic flow data [25] indicating a necessity to 
use some kind of viscous fluid description. A nonrelativistic viscous fluid is usually described by the first order 
Navier-Stokes equations. However, one needs their special relativistic generalization and this turns out to be 
acausal and unstable [26] (see also [27-29]). One therefore looks towards the extended, second order theories 
accepting all problems connected with their formulation and practical applications [30-36]. Physically, the 
difference is in that first order theories are based on the local equilibrium hypothesis, in which the independent 
variables are used, whereas in higher order theories the fluxes of the local equilibrium theory appear as 
independent variables. In particular, the entropy vector is quadratic in the fluxes, containing terms characterizing 
the deviation from local equilibrium. This situation plus our experience with the nonextensive formalism 
[10-12, 14-16] prompted us to investigate the simple nonextensive formulation of the perfect hydrodynamic 
model, a perfect q-hydrodynamics [37]. It turned out that this describes the experimental data fairly well. In 
addition, an apparently unexpected feature appeared, namely the possibility that (relatively simple and first 
order) perfect q-hydrodynamics can serve as a model of (second order and complicated in practical use) viscous 
d- hydrodynamics. This is the point we would like to discuss in more detail in this work. 

In the next Section we shall, for completeness, present the main points of q- hydrodynamics [37]. In Section 3 we 
propose a nonextensive/dissipative conjecture (NexDC), which allows us to connect ideal q-hydrodynamics with 
a d-hydrodynamics. Consequences of NexDC are discussed in Section 4 (entropy production) and in Section 5 
(transport coefficients). Section 6 contains the summary. 

2. Basic elements of q- hydrodynamics 

As in [37], we shall limit ourselves to a 1 + 1 dimensional baryon-free version of hydrodynamic flow. This is 
derived following Lavagno [38], in which a nonextensive version of the Boltzmann equation has been proposed 
and investigated. Because no external currents are assumed, this corresponds to a kind of perfect q-fluid. There 
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are two important points in [38]: (i) the Boltzmann equation is formulated for f q (x,p) distribution rather than 
for the usual f q (x,p); (ii) the usual molecular chaos hypothesis is now assumed in nonextensive form: 



hq\fq, fgi] = exp [In, f q + In, f ql ] where exp AX) = [1 + (1 - q)X] r=« , ln q (X) = 



(1) 



Here f q (x,p) is the g- version of the phase space distribution function, whereas h q [f q , f q i] is the q- version of the 
correlation function related to the presence of two particles in the same space-time position x but with different 
four-momenta p and pi, respectively. By postulating Eq. (1) we are, in fact, assuming that, instead of a strict 
(local) equilibrium, a kind of stationary state is being formed, which already includes some interactions and which 
is summarily characterized by a parameter q; very much in the spirit of [20-23] mentioned before. With this 
assumption we are already departing from the picture of the usual ideal fluid with its local thermal equilibrium, 
which is the prerogative of ideal hydrodynamics [26]. The consequences of this fact will be discussed below. 
The most important ingredient for further discussion is now the corresponding nonextensive entropy (g-entropy) 
current [37, 38]: 



d s p 



»{f q q^qfq- fq}- 



(2) 



It turns out that d^a q > at any space-time point, i.e., the relativistic local ff-theorem is valid in this case 
[37-39] cndnotc2 Assuming d M cr£ = one finds that (k B is Boltzmann constant): 



fq{x,p) 



» M <W ] 1/(1 " 9) 
1 q) k B T q (x)_ 



exp„ 



'k B T q {x) 



(3) 



where T q (x) is the g-temperature [37] and u%(x) is the g-hydrodynamic flow four- vector. Actually, one should be 
aware of the fact that there is still an ongoing discussion on the meaning of temperature in nonextensive systems. 
However, the small values of the parameter g — 1 deduced from data allow us to argue that, in the first approxima- 
tion, T q can be regarded as the hadronizing temperature in such a system (cf., [40] for a thorough discussion of the 
temperature of nonextensive systems). Finally, we get the q- version of the local energy- momentum conservation 
[38], 



d v T q uv (x) = Q, with T q ^(x) = 



exp 



k B T q (x 

In what follows we shall use covariant derivative notation in which the vector u M and tensor <? M " are defined as 



(4) 



«i„ = c/ M u + r Afl u and = d M gr + T£ M g + 



(5) 



by means of the Christoffel symbols, Y v Xll = ^g va (d^g^x + d\g al i — d^gxn). In this notation Eq. (4) reads: 
= + = [^PiKX-^WAT] =0, 



(6) 
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where A," = g^ v — u q u q . Here it was assumed that the g-modified energy-momentum tensor T q v can be 
decomposed in the usual way in terms of the g-modified energy density and pressure, e q = u ql /T q ^ v u qlJ and 
P q = —jjT^Aqitv, by using the g-modified flow u q (such that in the rest frame of the fluid u q — (1,0,0,0) and 
for q — y 1 it becomes the usual hydrodynamic flow u M ). Notice that Eq. (6) is formally identical to the perfect 
hydrodynamic equation but with all the usual ingredients replaced by their g-counterparts (perfect means here 
that there is nothing on the r.h.s. of Eq. (6)). In this sense we can speak about the perfect q-fluid mentioned before. 

Some remarks are in order before proceeding further. 

• When applied to multiparticle production processes, each hydrodynamic model is supplemented with three 
ingredients, which must be considered together with Eq. (6) (or its equivalent): (i) - initial conditions 
(IC) setting the initial energy density which is going to evolve hydrodynamically; (ii) - equation of state 
(EoS) reflecting the internal dynamics of the fluid considered and (m) freeze-out conditions describing 
transformation of the expanded fluid into observed hadrons. The problem is that all of them can, in 
principle, enter with their own intrinsic fluctuation pattern, i.e., with their own parameters q. In [37] where 
we provided preliminary comparison with experimental data, we have assumed, for simplicity, the same 
value of parameter q throughout the whole collision. It is still to be checked how good this assumption is. 
However, in our present discussion this point is unimportant. 

• Whereas in the usual perfect hydrodynamics (based on the BG statistics) entropy is conserved in the 
hydrodynamic evolution, both locally and globally, in the nonextensive approach it is only conserved locally. 
The total entropy of the whole expanding system is not conserved, because for any two volumes of the fluid, 
Vi, 2, one finds that Sq + S q V2 ' 7^ S q Vl<SV2 ^ (where S q V ^ are the corresponding total entropies) . This should 
be always remembered (albeit, strictly speaking, the hydrodynamic model requires only local, not global, 
entropy conservation). As a consequence of this fact, as we shall see below, contrary to the situation in 
usual perfect hydrodynamics, in the perfect g-hydrodynamics the entropy is produced (but not g-entropy). 

• To guarantee that hydrodynamics makes sense, there should exists some spacial scale L such that the 
volume L 3 contains enough particles. However, in the case when there are fluctuations and/or correlations 
characterized by some typical correlation length I for which we expect that I > L, one has to use nonextensive 
entropy S q and its (locally defined) density, s q (x) = S q /L . When formulating the corresponding q- 
hydrodynamics one takes the limit L — -> 0, in which case the explicit dependence on the scale L vanishes, 
whereas the correlation length I leaves its imprint as a parameter q. In this sense, perfect g-hydrodynamics 
can be considered as preserving causality and nonextensivity q is then related with the correlation length 
I (one can argue that, very roughly, q ~ l/L e g > 1, where L e g is some effective spacial scale of the q- 
hydrodynamics). If the correlation length I is compatible with the scale L e g, i.e., I w L e g, one recovers the 
condition of the usual local thermal equilibrium and in this case the g-hydrodynamics reduces to the usual 
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(BG) hydrodynamics. The above considerations were limited to the case of q > 1 only, for which, as was 
said before, the clear correspondence with fluctuations was found. In what concerns the case of q < 1, it 
seems that following [17], where q < 1 was found as a main factor closing the allowed phase space, we can 
at the moment only propose that it could probably correspond to the case where, for some reason, the scale 
L cannot vanish but must stop at some value I < L. In this case, analogously to what was said above, one 
could expect that, again, q ~ l/L e g, which this time would be smaller than unity. We shall not discuss this 
possibility further in this work. 

3. Nonextensive/dissipative conjecture (NexDC) 

As seen in Eq. (6), the structure of the perfect q-hydrodynamical flow is formally identical with the flow of ideal 
fluid described by the usual ideal hydrodynamics. We therefore regard the fluid described by Eq. (6) as a perfect 
q-fluid (in the sense already mentioned before). Our experience with applications of q statistics [12, 13] tells us 
that in cases of interest to us \q— 1| << 1. It is therefore tempting to simply expand the corresponding quantities 
in powers of \q — 1| and to only keep the linear term [23]. The result one gets looks promising, namely 

Tq V = ^ + (9 -IX" where (7) 

with T q =i being the usual energy-momentum tensor describing an ideal fluid in the BG approach and r q v repre- 
senting a viscous correction caused by the nonextensivity. However, in order for Eq. (7) to be valid in the whole 
phase space, the |1 — q\ (^?) 2 < 2 inequality must hold. This means that either such a procedure can be applied 
only to a limited domain of phase space, or that q = q(x,p), a possibility which is outside the scope of the present 
work. We must therefore proceed in a more general way. Let us formally decompose e q , P q and u q in Eq. (6) 
into the, respectively, extensive and nonextensive parts: 

e q (T q ) =e(T q ) + Ae q (T q ), P q (T q ) = P(T q ) + AP q (T q ), u%(x) = u»{x) + 5u^(x). (9) 

Actually this can be only done approximately because our extensive e and P still depend on T q , rather than on 
T. We tacitly assume that T q is not too far from T (so far they are independent parameters but later on we shall 
impose a condition on them, see Eq. (13) below). In Eq. (9) the four- velocity it M (x) is formally a solution of the 
equation to which Eq. (6) is transformed by using Eq. (9): 

\i(T q )uV - P{T q )A^ + 2W iM u ) + tt^] = 0. (10) 

One can see that it has form of a dissipative hydrodynamic equation [30-36] (here A M ' y = g^" — u^u" and 
A^B V) = \{A"B V + A U B^)) where e is the energy density, P is the pressure, W 1 denotes the energy or heat 
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flow vector, whereas 7r M " is the shear (symmetric and traceless) pressure tensor. They are defined as (using the 
angular bracket notation: a<^b v> % f [f (A£A£ + A£A$;) - |A M "A A(T ] a x b a ): 

e = s 9 + 3n, P = P q +n, ^ = ^ q [l+ 7 ]A^4, tt^ = + IIA^ = w q SufSuf, (11) 

and expressed in terms of IT (being a kind of a g-dependent bulk pressure), g-enthalpy w q and a g-dependent 
variable 7, 

II = ^w q [j 2 + 27], w q = e q + P q 7 = u M <5u£ = -^Suq^Su 1 ^. (12) 

Notice that when q — ► 1, the difference between the g and ideal hydrodynamic flows vanishes, <5£ — > 0, and with 
it also 7^0. It means that all dissipative fluxes of d-hydrodynamics (10) which are induced by the (/-flow, like 
W^, 7r M " and II, vanish in this limit as well, and one recovers the equations of the usual perfect hydrodynamics. 
The variable 7 is easier to handle and to calculate than the differences in flows (see Eq. (18) below for its explicit 
form). Notice that, whereas the time evolution of II is controlled by g-hydrodynamics (via the respective time 
dependencies of s q , P q and 7) its form is determined by the assumed constraints which must assure that the 
local entropy production is never negative (as is always assumed in the standard 2 nd order theory [31, 32] ). 

Now comes the crucial point. To finally link the usual g-hydrodynamics and its d counterpart, one has to fix 
somehow the temperature T q and the flow velocity field u q . This is done by assuming that there exists such a 
temperature T and velocity difference 5u^ that the following two relations are satisfied: 

P(T) = P q (T q ), e(T) = £q (T q )+m. (13) 

We call them the nonextensive/ dissipative relations, NexDC in short. Here e and P are the energy density and 
pressure as defined in the usual BG statistics, i.e., for q = 1. Eq. (13) provides a definite relation between T q 
and T, therefore, in what follows, we shall mainly use T for a description of dissipative effects (except of some 
expressions when it is easier to keep both T and T q , but always with the understanding that, because of Eq. (13) 
they are not independent). It is now straightforward to show [37] that in this case one can transform Eq. (10) 
into the equation of the usual d-hydrodynamics: 

{ £ (TKV - [P{T)+Il]A^ + 2W^u v) + tt""} =0. (14) 

This completes a demonstration of equivalence of the perfect g-hydrodynamics represented by Eq. (6) and its 
d-hydrodynamical counterpart represented by Eq. (14). It means therefore that the perfect g-fluid is nothing 
but a viscous fluid which satisfies the d-hydrodynamic equation Eq. (14). 
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With the bulk pressure II given by Eq. (12) and using the NexDC relations (13) one can express the q-enthalpy 
w q by the usual enthalpy, w = Ts = e + P, and the (/-dependent variable 7: 



= e q (T q ) + P q (T q ) = 



e(T) + P(T) 



[1 + 7 P [1 + T ]2 



(15) 



Notice that Eq. (13) leads to the following important relations between the heat flow vector W /M , the pressure 
tensor 7r'" / and the bulk pressure II: 



W M W„ = -3U.W, if"W v = -2W 



7Tu^7r M!/ = 611 



(16) 



As mentioned before, we expect that in all cases of interest to us \q — 1| << 1. This means then that we can 
regard the state characterized by f q (x,p) from Eq. (3) as some stationary state existing near equilibrium. This 
near equilibrium state is then defined by the correlation function h q in Eq. (1), for which the energy momentum 
tensor can be divided into two parts: the energy-momentum tensor of the usual ideal fluid, T q Z\ = TJtf , and the 
g-dependent remaining <5T M " (the meaning of all components is the same as in Eqs. (11) and (12) and ST^ V — > 
when q — > 1): 

Tf = (e q + P 9 )« - P q g" v d = T£ + 5T" V where 5T» V = -IIA^ + 2W ( "/» + vr"" . (17) 
Using now Eq. (13) we get from Eq. (15) that 



7 = y/l + Se q - 1 



where 5e„ 



e(T) - e q {T q ) 
£ q {T q ) + P q (T q ) 



(18) 



This relation allows us to connect the velocity field u q (which is solution of the g-hydrodynamics given by Eq. (6)) 
with the velocity field u (which is solution of the dissipative hydrodynamics given by Eq. (14)). To make it more 
transparent, let us parameterize these velocity fields by using the respective fluid rapidities a q and a (the metric 
used is g» v = (1,-1/t 2 )): 



u q (x) 



cosh (a q — t?) , — sinh (a q — 77) 

T 



cosh (q — 77) , — sinh (a — 77) 
r 



(19) 



Because 7 = 7t M 5u^ = cosh(a 9 — a) — 1 the connection between u and u q is given by 



sh(a 9 — a) = y/l+~6e 



a = a„ — lot 



; (e q + sjl + Se q ^ 



(20) 



(we abandon the other solution, a = ct q + log i^e q + %/T-j- Se q ^j , because it leads to the reduction of entropy, not 
to its production, i.e., for it [su^]-^ < 0, for q > 1). 
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4. Entropy production in (/-hydrodynamics 

One of the important implications of the NexDC conjecture is that in the ideal g-hydrodynamics one observes 
entropy production. Taking the covariant derivative of Eq. (17) and multiplying it by u v one gets 

u v T^ = T[su^ + u„5T» v = =► = -^ST.^. (21) 

This means that, although in the ideal g-hydrodynamics the g-entropy is conserved, i.e., [sqU^]-^ = 0, the usual 
entropy is produced; the ideal g-fluid is therefore a kind of usual viscous fluid. This is illustrated in Fig. 1 where the 




Figure 1 . Illustration of different aspects of the entropy production, [su^] ;M > 0. Two left panels show production of entropy taking 
place in different parts of the longitudinal phase space characterized by the rapidity variable, y — (1/2) ln[(i + z)/(f. — z)] , 
observed at different evolution moments r. First panel shows results for the initial conditions given by fixed maximal 
initial energy density e^ in \ second panel shows results for fixed width of the initial energy distribution (assumed in 
gaussian form) a - sec text for other details (cf. , also [37]). Right (third) panel shows how entropy is produced in the 
central region (for y — 0, first type of initial conditions were used here) for different values of the parameter q and for 
different evolution moments r. This entropy production results in the corresponding growth with q of the multiplicity 
of produced secondaries, N = 11 • (7 + 10q) (cf. , Fig. 8 in [37]). 

expected entropy production, as given by Eq. (21), is shown. All curves presented in Fig. 1 were calculated using 
the g-hydrodynamical evolution described in [37] for q — 1.08 with the g-dependent initial conditions (q-lC). They 
were given assuming a g-gaussian shape (in rapidity) of the initial energy density out of which the hydrodynamic 
evolution started. Two types of (g-IC) were considered: in the first the maximal initial energy density was fixed 
as e' 1 ™' = 22.3 GeV/fm 3 (forcing the width of the g-gaussian to be equal to a — 1.28), in the second the width 
of the g-gaussian was fixed as a = 1.25 (forcing the initial energy density to be equal to e' ln ' = 27.8 GeV/fm 3 ). 
Because both (g-IC) give reasonable fits to experimental data, they introduce only very small differences in the 
entropy production, cf., the left part of Fig. 1. These (g-IC) were accompanied by the g-dependent equation of 
state (g-EoS) for the relativistic pion gas. Finally, the usual Cooper-Frye prescription of the freeze-out (but with 
the corresponding distribution being the g-exponent rather than the usual one [37]) was used. All calculations 
were performed for Au + Au collisions at ^/snn = 200 GeV) (see [37] where it was shown that g-hydrodynamics 
with such g and with the respective g-IC and g-EoS reproduces dN/dy and px distributions observed in RHIC 
experiment for Au+Au collisions at ^/snn = 200 GeV energy). Notice that suj^ > for large n region at any r 
(but especially for the early stage of the hydrodynamical evolution). It supports therefore a dissipative character 
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of the g-hydrodynamics mentioned before and leads us to the conclusion that the equilibrium state generated in 
heavy-ion collisions may, in fact, be the g-equilibrium state, i.e., some stationary state near the usual equilibrium 
state already containing some dissipative phenomena. Notice that the total multiplicity, which is usually treated 
as a measure of entropy, increases with q as expected, namely N = 11 • (7 + 10g). 

5. Transport coefficients of q- hydro dynamics. 

The other implication of the NexDC conjecture is that there are some transport phenomena in an ideal g-fluid. 
We shall now identify the corresponding transport coefficients and compare them with the similar coefficients in 
the usual viscous d- hydrodynamics. Out of a number of different formulations of d-hydrodynamics [27-36] we 
shall choose the 2nd order theory of dissipative fluids presented in [31, 32]. It does not violate causality (at least 
not the global one over a distant scale given by the relaxation time) and it contains some dissipative fluxes like 
heat conductivity and bulk and shear viscosities. We shall now see, to what extent dissipative fluxes resulting 
from our g-hydrodynamics can be identified with the heat conductivity and with the bulk and shear viscosities 
introduced in [31, 32], and what are the resulting transport coefficients. 

Let us start with the most general form of the off-equilibrium four-entropy current cr M which in our baryon-free 
scenario takes the following form [31, 32]: 

a" = P{T)P" + (3„T» V + Q^iST"") with /3 M = ^, (22) 

where ST' 11 ' is defined in Eq. (17). The function Q M = Q M (ST 1111 ) characterizes the off-equilibrium state, which 
in our case is induced by the nonextensive effects and therefore depends on g, Q — > when q 1. Using the 
NexDC conjecture, Eqs. (13), and Eq. (15) our g-entropy current is 

<(ee <) = su» + ?p + Q» where Q£ = Q^ x = X + ^ J . (23) 

Our near equilibrium state is thus characterizes by the function in which 

T q V w 1 + 7 T q 

(temperatures T and T q are not independent but connected, as was stated before, by the NexDC conjecture 
(13)). We stress that it is given by the dissipative part of our g-system, the one that leads to the increase of the 
usual entropy, cf. Fig. 1 (the g-entropy, as was discussed before, remains, however, strictly conserved, a q .^ — 0). 

The most general algebraic form of for d-hydrodynamic that includes dissipative fluxes up to second order is 
[32]: 



[-ATI 2 + ^W V W V - f}2*„x* vX ] „ aoW . a^W v 
^2nd — 2T u T T ' " ' 
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where /3j=o,i,2 are the corresponding thermodynamic coefficients for the, respectiveiy, scaiar, vector and tensor 
dissipative contributions to the entropy current, whereas Oj=o,i are the corresponding viscous/heat coupiing 
coefficients. The corresponding expression caicuiated using the NexDC conjecture is: 



= r 2n d SU M + Tit 



' T 



where 



3^i n _ (A) + 6/3 2 ) n2 
2 2w 



and 



T, 



(a + 2 ai )n. 



(26) 
(27) 



This is given by the second order polynomiai in the buik pressure II. Therefore, it is naturai to expect (and this 
wiil be our assumption) that the most generai entropy current in the NexDC approach is: 



r(n)su" + T(n) 



T 



(x + + ( X - 



where 



x = (r + T)/2 



and 



T ' 

(r - T)/2 



(28) 
(29) 



and where T, T are (in general infinite) series in powers of the bulk pressure II. In this sense Qf ull can be regarded 
as the full order dissipative current in (/-hydrodynamics. In general one has entropy production/reduction, i.e., 
crtt ^ 0, however, in the case when r(II) = T(II) = \ one has cr£ ;jJ = 0. Out of the two possible solutions for 
(r, T) only one is acceptable, 



2? 



3n 

w 



2 7 T 
1+7^ 



T 



(30) 



because only for it it/jQ^n < (i.e., the entropy is maximal in the equilibrium [32], this is because (T — T q )/T q is 
always positive for q > 1 [37]). In this way we finally arrive at the following expression for the full order dissipative 
entropy current emerging from the NexDC approach endnote3 : 



su + 



w 



2T 

Ta 



1 - \ 1 ~ 



311 



su + 



2(T 



(31) 



Limiting ourselves to situations when T/T q « 1 and neglecting terms higher than C(3II/ui) 2 , one obtains that 



QiuU 



m 

w 



4 I w 



(32) 



Comparing now Eqs. (27) and (32) one gets that /3i 



Po + 6/3 2 



and Qfo + 2a i = 0. Since in the 



Israel-Stewart theory [31] the relaxation time r is proportional to the thermodynamic coefficients /3o,i,2, it is 
naturally to assume that in the case of NexDC it is proportional to the inverse of the enthalpy, r oc l/w (notice 
that for the classical Boltzmann gas of massless particles /3 2 = 3/w [32, 36]). 
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Let us now see what kind of bulk and shear viscosities emerge from the NexDC approach. To this end let us write 
the full order entropy current Eq. (28) in the following form: 



[(1 + x)$H;m + K* M Lm wh ere = si/ 1 



T 



and = su" - ~. (33) 



Because conservation of the g-entropy, a^.^ = 0, is equivalent to [(1 + x) < J >m ];m = 0> therefore, using Eq. (16) one 
gets that ^ = ~m^u" + ^gsTr"", therefore 



where : 



X u = --, 



n 



n 



n t £ 



Y a = 



n 



2 1 1 



YyW v 

2n 



(34) 
(35) 

(36) 
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Eq. (16) allows to eliminate the term proportional to the heat flow, 



Finally one obtains 



'full;H 



I i TT^ - 



n 2 

CT 



2rjT 



>0, 



(37) 



where we have introduced the usual bulk and shear viscosities, £ and r\. Notice that, because of Eq. (16) one 
avoids the explicit contribution to the entropy production coming from the heat flow, which is present in 
Eq. (34) when one discusses a baryon free fluid, in which case the necessity to use the Landau frame would appear. 
As one can see, Eq. (33) is covariant and therefore it does not depend on the frame used. One arrives at our main 
result: the sum rule connecting bulk and shear viscosity coefficients (expressed as their ratios over the entropy 
density s), see Fig. 2: 



1 3_ 

C/ s ri/s 



"full;, 
IP 



(38) 



This is as much as one can get from the q-hydrodynamics alone. To disentangle this sum rule, one has to add 
some additional input. Suppose, therefore, that we are interested in an extremal case, when a total entropy is 
generated by action of the shear viscosity only. In this case one can rewrite the first part of Eq. (37) as 



[" en {wuXjCx) + ~ Z A ' 



(39) 



and arrive at 



V = 7(7 + 2) 
s (7+ I) 2 



II 



— SU A a 



(40) 



The predictions of Eq. (40) are shown in Fig. 3 for different values of the parameter q and for different rapidities. 
They are also confronted with the known result on ri/s, provided by the popular Ads/CFT conjecture [41], that 
v/ s > 1/47T. Assuming the validity of this limitation we can use Eq. (38) only in the region where the r.h.s. 
of Eq. (40) is smaller than (or equal to) 1/47T (in which case we put rj/s = l/47r), otherwise (because of our 
assumption that the total entropy is generated by the shear viscosity only) we have to put £/s = and use 
Eq. (40) to evaluate 77/ s. The results for £/s and rj/s are shown in Fig. 3. Notice that, when the r.h.s of Eq. (40) 
approaches 1/47T, £/s given by Eq. (38) approaches infinity. To avoid such a situation, rj/s should start to increase 
at higher temperatures, for example at T > 75 MeV cndnotc4 . 



6. Summary 

To summarize, we have discussed dissipative hydrodynamics from a novel point of view. This is provided by 
the nonextensive formulation of the usual perfect hydrodynamical model recently proposed by us [37]. Such a 
model can be solved exactly and contains terms which can be interpreted as due to some dissipative effects. 
They can be identified and are expressed by the nonextensivity parameter q of the Tsallis formalism applied 



Noncxtensive perfect hydrodynamics - a model of dissipative relativistic hydrodynamics? 
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Figure 3. Upper panels: the ratio of the shear viscosity over the entropy density, r]/s, as function of temperature T, calculated 
using Eq. (40) for different values of q at the mid-rapidity region y — (two upper-left panels) and for rapidity y — 3 
(two upper-right panels). The same version of hydrodynamical model was used as in Fig. 1 with the first type of 
initial conditions. Notice that for q — > 1 this ratio vanishes, as expected, proving therefore correctness of our numerical 
calculations. The lower two panels summarize the rapidity dependence of both the shear and the bulk viscosity ratios, 
r]/s (left panel) and £/s (right panel). The rj/s is calculated using Eq. (40), i.e., assuming Q/s — 0, and presented only 
above the limit value rj/s — 1/4-7T found in the AdS/CFT approach [41]. The right panel shows the bulk viscosity, C/s, 
calculated in this limit, i.e., assuming rj/s — 1/Air and using Eq. (38). 



here. This finding was used to propose a possible full order expression for the dissipative entropy current a^ ull 
resulting from the nonextensive approach. The corresponding bulk and shear transport coefficients resulting 
from q-hydrodynamics are connected by a kind of sum rule, Eq. (38). They were calculated for some specific 
simplified case, cf. Fig. 3. We close by noticing that there is still some uncertainty in the equation of state used 
in our numerical example (for example, in what concerns the role of the possible QGP phase transition, which, 
when included in EoS, could considerably affect rj/s presented in Fig. 3). We plan to address this subject elsewhere. 

From our point of view it is interesting to observe that (at least some) remedies proposed to improve the formu- 
lation of d-hydrodynamic [26, 29] (like the use of some induced memory effects) introduce conditions which in 
statistical physics lead in a natural way to its nonextensive version described by the parameter q [9]. It is then 
natural to expect that a nonextensive version of the hydrodynamical model with the nonextensivity parameter 
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q could help us to circumvent (at least to some extent) the problems mentioned. It is because equations of the 
perfect nonextensive hydrodynamics (or perfect g-hydrodynamics) can be formally solved in an analogous way as 
equations of the usual perfect hydrodynamics [37]. However, it turns out that, from the point of view of the usual 
(extensive) approach the new equations contain terms which can be formally identified with terms appearing in 
the usual dissipative hydrodynamics (d- hydrodynamics). Although this does not fully solve the problems of d- 
hydrodynamics, nevertheless it allows us to extend the usual perfect fluid approach (using only one new parameter 
q) well behind its usual limits, namely toward the regions reserved so far only for the dissipative approach. 
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Notes 

[42] One should keep in mind, however, that there are dualities in the non-extensive approach, i.e., that both 
q and 1/q can be used as the nonextensivity parameter depending on the normalization of original or 
g-powered probabilities. Also, when considering the particle-hole symmetry in the g-Fermi distribution, 
f(—E, T, (x, q) = 1 — f(E, T, — /i, 2 — q), in a plasma containing both particles and antiparticles both q and 
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2 — q occur (/i denotes the chemical potential here). These dual possibilities warn us that a theory requiring 
that only q > 1 has physical meaning is still incomplete [9]. These points deserve further considerations 
which are, however, going outside the scope of this paper. On the other hand, in the present paper we 
are not considering plasmas containing both particles and antiparticles but only massive pions which are 
assumed to obey the g-Boltzmann distribution. 

[43] Actually, this form of q-entropy current differs slightly from that in [38]. The reason is that only with 
such form we can, at the same time, both satisfy the _ff-theorem and reproduce thermodynamical relations 
resulting in our (/-enthalpy, cf. Eq. (15), ((cf. also Eqs. (17a) and (17b) in [37]), which are crucial in all 
further derivations and which were not addressed in [38]). 

[44] One must keep in mind therefore that, although to define the entropy current in q-statistics we require 
only that q > 0, the NexDC correspondence requires in addition that q > 1 to be consistent with d- 
hydrodynamics [32]. Only then, as witnessed by Fig. 1, constant (nonzero) initial Tsallis entropy results 
in increasing BG entropy (demonstrating itself in the total multiplicity increasing with (q — 1)). For q < 1 
we could have decreasing BG entropy from some positive initial value (and, accordingly, total multiplicity 
decreasing with (q — 1)). This point is, however, not totally clear at present and we plan to address it 
elsewhere. 

[45] It should be kept in mind that we have so far obtained only the relations between the Israel-Steward 
coefficients (/3o,i,2 and ao,i) an< A n °t their individual values, therefore we cannot compare our results to 
those of the ideal Boltzmann gas case. This is due to the tensor form of relations (16) and to the fact 
that the original perfect g-hydrodynamics (Eq. (4)) does not contain any natural space-time scale. It is 
then natural that NexDC conjecture does not introduce per se any definite relaxation time or viscous-heat 
coupling length scale. Also, in examples shown here, only EoS for pionic gas was used as in [37]. Although 
it was shown there that such EoS depends only very weakly on the parameter q, it remains to be checked 
whether this is also true for a more realistic EoS with quarks and gluons and in the vivinity of the QGP 
— » hadronic matter phase transition. We plan to address this point elsewhere. 



